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Using fully kinetic simulations, we study the scaling of the inflow speed of collisionless magnetic 
reconnection in electron-positron plasmas from the non-relativistic to ultra-relativistic limit. In the 
anti-parallel configuration, the inflow speed increases with the upstream magnetization parameter 
a and approaches the speed of light when a > 0(100), leading to an enhanced reconnection rate. 

In all regimes, the divergence of the pressure tensor is the dominant term responsible for breaking 
the frozen-in condition at the x-line. The observed scaling agrees well with a simple model that 
accounts for the Lorentz contraction of the plasma passing through the diffusion region. The results 
demonstrate that the aspect ratio of the diffusion region, modified by the compression factor of 
proper density, remains ~ 0.1 in both the non-relativistic and relativistic limits. 
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Introduction- Magnetic reconnection is a process that 
changes the topology of magnetic fields and often leads 
to an explosive release of magnetic energy in nature. It is 
thought to play a key role in many energetic phenomena 
in space, laboratory and astrophysical plasmas [1]. In 
recent years, relativistic reconnection has attracted in- 
creased attention for its potential of dissipating the mag- 
netic energy and producing high-energy cosmic rays and 
emissions in magnetically dominated astrophysical sys- 
tems [2], such as pulsar winds [3-5], gamma-ray bursts 
[6-8] and jets from active galactic nuclei [9-11]. However, 
many of the key properties of magnetic reconnection in 
the relativistic regime remain poorly understood. While 
early work found the rate of relativistic magnetic recon- 
nection may increase compared to the nonrelativistic case 
due to the enhanced inflow arising from the Lorentz con- 
traction of plasma passing through the diffusion region 
[12, 13], it was later pointed out that within a steady- 
state Sweet-Parker model [14, 15] the thermal pressure 
within the current sheet will constrain the outflow to 
mildly relativistic conditions where the Lorentz contrac- 
tion is negligible [16], and a relativistic inflow is therefore 
impossible. Recently, the role of temperature anisotropy 
[17], inflow plasma pressure [18], two-fluid [18], inertia 
effects [19] and mass ratio [20] have been considered. All 
existing theories are generalizations of the steady-state 
Sweet-Parker or Petschek-type [21] models, which do not 
account for the mechanism that localizes the diffusion re- 
gion and determines the reconnection rate in collisionless 
plasmas. Meanwhile, a range of reconnection rates are 
reported in computational works with different simula- 
tion models and normalization definitions [18, 20, 22-25]. 
However, the scaling of the rate has yet to be determined 
and the kinetic physics of the diffusion region is poorly 
understood in the relativistic limit. 

In this work, a series of two-dimensional (2D) full 
particle-in-cell (PIC) simulations have been performed to 
understand the properties of reconnection in the magneti- 


cally dominated regime. It has been argued that electron- 
positron pairs are relevant in highly energetic astrophys- 
ical environment, such as pulsar winds [4, 26] and extra- 
galactic jets [27], hence in this Letter we limit our study 
to mass ratio mi/m e = 1. The magnetization parameter 
can be defined as the ratio of the magnetic energy den- 
sity to the plasma energy density, a = B 2 /(Attic) with 
enthalpy w = 2 n'mc 2 + [T / (T — 1)]P' . Here T is the ratio 
of specific heats and P' = n’ (T' e + T p ) the plasma ther- 
mal pressure in the rest frame. The shear Alfven speed 
is Va = c[cr / (1 + cr)] 1 / 2 [18, 28-30]. In this Letter, the 
primed quantities are measured in the fluid rest (proper) 
frame, while the unprimed quantities are measured in the 
simulation frame unless otherwise specified. As pointed 
out in Ref. [16], if a simple pressure balance P' ~ B 2 /8tt 
is satisfied across a steady- state Sweet-Parker layer, this 
will constrain the effective a ~ 0(1), and thus restrict 
the inflow speed V{ n <C c. However, we demonstrate the 
development of relativistic inflows when the upstream 
a > 0(100) (for the first time) in fully kinetic simula- 
tions. A simple model based on the underlying idea of 
Blackman and Field [12] is presented to explain the scal- 
ing of the inflow speed and normalized reconnection rate. 
It is well known that the normalized collisionless recon- 
nection rate, R = Vinl^Ax ~ 0.1, in the non-relativistic 
limit can be estimated by the aspect ratio of the diffusion 
region, but the precise physics that determines this value 
remains mysterious [31]. Here Va x is the Alfven wave 
velocity in the outflow direction. Interestingly, the simu- 
lation results in this study suggest that this aspect ratio 
(modified by the compression factor of proper density at 
the inflow and outflow) of ~ 0.1 persists in the relativistic 
regime. In addition, we analyze the relativistic general- 
ization of Ohm’s law [32], and identify the importance of 
the pressure tensor and the time- derivative of the inertial 
term in breaking the frozen-in condition. 

Simulation setup- The majority of simulations in this 
letter start from a relativistic Harris sheet [22, 33-35]. 
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The initial magnetic field B = F> 0 tanh(z/A)x + B g y cor- 
responds to a layer of half-thickness A with a shear an- 
gle 0 = 2ta,n~ 1 (Bo/B g ). Each species has a distribution 
fh oc sech 2 (z/ X)ex.p[—^d{lLmc 2 -\-mVdUy) /T'} in the sim- 
ulation frame, which is a component with a peak density 
n' 0 and temperature T' boosted by a drift velocity ±Vd in 
the y-direction for positrons and electrons, respectively. 
Here u = is the spacelike components of 4- velocity, 
7 l = 1/[1 — (^/c) 2 ] 1 / 2 is the Lorentz factor of a par- 
ticle, and 7 d = 1/ [1 — (V^/c) 2 ] 1 / 2 . The drift velocity 
is determined by Ampere’s law ci?o/(47rA) = 2e7^noV^. 
The temperature is determined by the pressure balance 
E>o/(87r) = 2n' 0 T' . The resulting density in the simula- 
tion frame is no = Jan'o- In addition, a non-drifting back- 
ground component fb oc exp (—7 ^me? fTf) with a uniform 
density n& is included. The simulations are performed us- 
ing VPIC [36] , which solves the fully relativistic dynamics 
of particles and electromagnetic fields. Densities are nor- 
malized by the initial background density n&, time is nor- 
malized by the plasma frequency uj pe = (47rn5e 2 /m e ) 1 / 2 , 
velocities are normalized by the light speed c, and spatial 
scales are normalized by the inertia length d e = c/uj pe . 
Although commonly used, the relativistic Harris sheet 
may not be generic. To test the sensitivity of our re- 
sults to the initial sheet equilibrium, a force-free config- 
uration [25, 37] was also included, with magnetic profile 
B = Hotanh(T/A)x + [B g + Hosech(^/A)] 1 / 2 y, and uni- 
form density and temperature. Particles in the central 
sheet have a net drift V p = — V e to satisfy Ampere’s 
law. All simulations use 100 — 200 particles per cell for 
each species (supplement D shows a convergence study). 
The boundary conditions are periodic in the x-direction, 
while in the z-direction the boundary conditions are con- 
ducting for fields and reflecting for particles. A local- 
ized perturbation with amplitude B z = 0.03Ho is used 
to induce a dominant x-line near the center of simula- 
tion domain. The simulation parameters for the var- 
ious runs considered in this Letter are summarized in 
Table 1. Our primary focus in the following section 
is the case Harris-4 which illustrates the dynamics in 
the transition to the limit with relativistic inflows (i.e., 
Vi n ~ c). The domain size is L x x L z = 384 d e x 384d e 
with 3072 x 6144 cells. The half-thickness of the initial 
sheet is A = d e , = n' 0 , Tb/mc 2 = 0.5, B g = 0 and 
uj pe /Vt C e = 0.05 where D ce = e5 0 /(m e c) is the cyclotron 
frequency. The upstream magnetization parameter based 
on the reconnecting component is a x = Bq/ (4ttw) = 
(D ce /u; pe ) 2 /{ 2[1 + (T/(T - l))(T 6 /mc 2 )]}, which is 88.9 
with T = 5/3. For cases with Tb/mc 2 > 1 in Table 1, we 
use T = 4/3 [38, 39]. 

Simulation results- Fig. 1(a) shows the structure of 
current sheet in the nonlinear stage, where the current 
density concentrates within a layer with a half-thickness 
^ d e . This thickness appears to be independent of the 
initial sheet thickness, and scales with the inertial length 
based on the asymptotic background density (rib). As 



FIG. 1: Results from case Harris-4 (a = a z = 88.9) during the 
fully nonlinear phase at 566.4/cj pe showing (a) the electron 
out-of-plane speed V ey , (b) outflow speed V ex with cut at 
z = 0, (c) inflow speed V ez with cut at z = —3.5 d e , (d) a 
closeup showing the non-ideal electric field E y + (V e x B)„ 
inside the green-dashed box depicted in (c). The non-ideal 
electric field is positive in between the horizontal white curves. 
Black contours are flux surfaces in (a)-(d). 


shown in Fig. 1(b), the outflow velocity approaches ~ c, 
while in Fig. 1(c) the peak inflow speed is ^ 0.65c. 
Note that these relativistic inflows also penetrate deeply 
across the magnetic separatrix into the d e -scale sheet 
in the downstream region (\x\ ~ 50 d e ). In addition, 
the simulation shows a rapid growth of secondary tear- 
ing modes, not only around the major x-line, but also 
along the concentrated current sheet that extends into 
the outflow exhausts. Interestingly, the secondary tear- 
ing mode appears considerably shorter spatially in com- 
parison with those in the non-relativistic regime. As 
shown in the blow-up (Fig. 1(d)), a magnetic island 
at (x ~ —2 d e ,z = 0) is immersed inside the region 
where the frozen-in condition is broken, and it has a 
size ^ 3 d e x 2 d e , implying that the secondary tearing 
mode grows for wave vectors k x S > 1. Here 5 is the 
half-t hickness of the intense nonlinear current layer. In 
contrast, the initial tearing mode based on the relativis- 
tic Harris equilibrium is still constrained by k x X < 1 
(i.e., from relativistic energy principle) [40], as in the 
non-relativistic limit. A temperature anisotropy [41, 42] 
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or the velocity shear associated with the outflow jet [43- 
45] may change the stability criterion, however, to resolve 
this issue in the relativistic regime is beyond the scope 
of this Letter. Fig. 1(d) shows that the non-ideal electric 
field is also concentrated in a region \z\ < d e . However, 
the frozen-in condition starts to fail inside a wider layer in 
between the horizontal white curves, which may be due to 
a larger effective inertial scale based on a smaller density 
at \z\ > d e (see the density cut in Fig. 2(a)). Fig. 2(a) 
shows that the inflow velocity V ez reaches its maximum 
~ 0.65c at the location where frozen-in starts to fail (i.e., 
marked by the green circle on the E y + ( V e x B ) y curve). 
The profile of V ez is rather flat in between this location 
and z = d e . Motivated by this observation, we use the lo- 
cal magnetic field B x , u at this location (z ~ 3.5 d e ) to nor- 
malize the reconnection electric field E y , and the normal- 
ized electric field traces the evolution of the peak inflow 
velocity well (Fig. 2(b)), as expected. Since in this rela- 
tivistic regime Va x ~ c, these two quantities are equiv- 
alent measurements of the normalized reconnection rate 
as discussed in the following section. The original peak 
density at the center of the sheet is no = 7^ + 1 « 11. 
This peak density drops significantly from 11 to ~ 2 and 
the density along the symmetry line (z = 0) remains 
~ 2 — 4, except inside secondary islands. The density 
ratio between the region immediately upstream to the 
x-line is ~ 2. 5/0. 3 = 8.3. These numbers will be used 
to estimate the compression factor in the following sec- 
tion. Per Ampere’s law, the density changes inside this 
d e -scale layer require a reduction of the local magnetic 
field since the motion of the current carrier is limited by 
the speed of light [46] . 

To examine the mechanism of flux breaking, we em- 
ploy the relativistic generalization of Ohm’s law E + 
V e x B + (l/en e ) V • P e + (m e /e)(d t U e + V e • VU e ) = 0. 
Here U = (1/n) f d 3 uuf, and the fluid velocity is V = 
(1/n) f d 3 uwf. The pressure tensor, P = f d 3 uvuf — 
nVU, defined in this manner reduces to the standard 
definition in the non-relativistic regime [32]. Each term 
along the vertical cut in Fig. 1(d) is plotted in Fig. 2(c). 
There are strong oscillations in both V -P e and V e • VU e , 
which largely cancel each other. In comparison, the mag- 
nitude of the non-ideal electric field E y + (V e x B) y is 
much smaller. In Fig. 2(d), we examine the region around 
the neutral point, which demonstrates that V e - VU e van- 
ishes at z = 0 since the neutral point coincides with the 
stagnation point in this symmetric configuration. The 
thermal pressure term, V • P e , balances the non-ideal 
electric field at the x-line while the time- derivative of 
the inertia d t U e remains small at this time [32], consis- 
tent with the study in the non-relativistic limit [47-50]. 
However, the intense current layer is strongly unstable to 
secondary tearing modes, similar to the non-relativistic 
limit [51]. The time-derivative of inertia df\J e becomes 
finite positive when the d e -scale current layer extends in 



FIG. 2: In (a), n e , E y + (V e xB)^ and \ V ez \ along the vertical 
cut shown in Fig. 1(d). \V ez \ is scaled by the right axis. The 
green circle marks the location where the frozen-in condition 
starts to fail; In (b), the evolution of the normalized reconnec- 
tion electric field E y /B XjU and the peak V ez near the major 
x-line at min (Ay) along z = 0. Here A y is the y-component 
of vector potential; In (c), quantities of Ohm’s law along the 
vertical cut shown in Fig. 1(d); In (d), the blow-up of (c) near 
the magnetic neutral point. 


length, and 9tU e becomes finite negative (i.e., contribut- 
ing to reconnection) when a secondary tearing starts to 
emerge in a sufficiently long layer (supplement B). 

Simple model- While previous theories [12, 13, 16, 19] 
generalize the Sweet-Parker [14, 15] or Petschek’s [21] 
models into the relativistic regime, here we simply an- 
alyze the conservation of mass including the influence 
of the Lorentz contraction over a control-volume of size 
L x 5, 


— hout'7 out^ out^ • (1) 

Here the subscript “m” and “ ouf indicate the inflowing 
and outflowing plasmas, respectively. Given a magnetic 
shear angle, the outflow is limited by the upstream Alfven 
wave velocity projected into the x-direction [31], 

V out = V Ax = cJEE. (2) 

V 1 + o' 

Here the upstream magnetization parameter is a = 
<j x + cr g with (j g = Bg/ (8 ttw) accounting for the contri- 
bution from the guide field. The effective Lorentz factor 
based on the bulk flows is 7 out = 1/ [1 — (V ou t/ c) 2 ] 1 / 2 = 

[(1 + <t)/( 1 + <7g)] 1/2 and 7 in = {(1 + a) /[l + a - 

a x (V in /V A ,) 2 ]} 1/2 . 

Working through the algebra, the peak inflow velocity 
can be determined with only one free parameter, r n '5/L , 

Vin = ( / Ox 

c \ n LJ Y 1 + <jg + (r n 'S/L) 2 a x ’ 
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FIG. 3: Scaling of the inflow velocity Vi n /c and the normal- 
ized reconnection rate R as a function of a x for cases with 
B g — 0 in the left; as a function of B g /Bo for cases with 
<j x — 88.9 in the right. Diamonds are measurements of runs 
in Table 1, green-dashed curves are predictions based on dif- 
ferent value of r n /5/L as marked on the plots. 


where r n > = n' out /n' in is the proper density ratio of the 
outflow to inflow. The compression factor is 


'Bout / lT<7 

n in r ” y 1 + Vg + (,r n >5/L) 2 <j x ’ 


and the normalized reconnection rate is 


( 4 ) 


„ _ Vin _ f 5_\ I 1 + <7 

_ Vax \ L ) y 1 + (Tg + {r n :5/L) 2 (j x ' 


( 5 ) 


Note that R differs from Vi n /c by a factor c/Va x • From 
the frozen-in condition, the normalized rate can also be 
written as R = (c/Vax)E v / B x , u - In the limit of Va x G 
then R ~ Vin/c ~ E y /B x , u as shown in Fig. 2(b). 

With the assumption of r n 'S/L = 0.1, as in the non- 
relativistic limit, Eqs. (3) and (5) immediately give 
R ~ Vi n /c = 0.69, consistent with the observed values for 
the case discussed. By comparing the measured compres- 
sion factor ~ 8.3 in Fig. 2(a) and n ou t/ui n = 6.9r n / from 
Eq (4), this implies that r n > ~ 0(1) and therefore the as- 
pect ratio 8/L ~ 0(0.1). The aspect ratio of the intense 
Ey + (V e x B) y layer shown in Fig. 2(d) seems to be con- 
sistent with this idea, however, a definite measurement 
is difficult because of time-dependency (more in supple- 
ment C). To further test these predictions, a series of runs 
were performed over a wide range of parameters (listed 
in Table 1). The measurement of Vi n /c and R are shown 
in Fig. 3 as diamonds, which agree closely with the pre- 
dicted scaling based on r n '8/L = 0.1. This suggests that 
the modified aspect ratio of the diffusion region persists 
during the transition from the non-relativistic to strongly 
relativistic regime. With a larger cr x , both the outflow 
and inflow speeds become closer to the speed of light. 
For anti-parallel initial conditions {i.e., a g = 0), both 
Vi n /c and R approach unity only when a x > 0(100), 


as shown in Fig. 3 (a)-(b), a condition obtained by de- 
manding (r n f 8 / L) 2 (i x = O.Olcr^ 1 in the denominator 
of Eq. (3) and (5). On the other hand, with a guide field 
Bg/Bo > 0(1), the outflow speed (Eq.(2)) becomes non- 
relativistic, the Lorentz contraction becomes negligible 
and the reconnection rate therefore goes back to ~ 0.1 
as shown in Fig. 3(d). To test the dependence on the 
choice of initial conditions, we have also performed an 
additional series of simulations using a force-free current 
sheet for initial condition [25]. The final states are simi- 
lar to those of initial Harris sheets and the measurement 
shown as blue diamonds in Fig. 3(a)- (b) follow the same 
trend, which demonstrates that the scaling in the nonlin- 
ear stage is determined solely by the upstream condition. 
Our model appears to explain the scaling of the normal- 
ized rate observed in two- fluid simulations of Zenitani et. 
al. [18] as well. Unfortunately, due to the complexity of 
evolution in the nonlinear phase, we are not able to pre- 
dict the reconnection rate normalized by the far upstream 
reconnecting component, which approaches a maximum 
of ^ 0.3 for cases with a x ~ 0(500) in the present study. 

Discussion- During the initial evolution of Harris-type 
current sheets, the pressure balance argument proposed 
in Ref. [16] restricts the inflow to Vi n <C c. However, 
at later times there are a variety of features which may 
break this argument. Firstly, the repeated formation of 
secondary plasmoids makes the diffusion region highly 
time-dependent. Secondly, the current density inside 
each of these plasmoids is much stronger than the cur- 
rent density within the diffusion regions between the plas- 
moids. This redistribution of current alters the structure 
of the reconnection layer and leads to strong variations 
in the reconnecting component of the upstream magnetic 
field. As a result, the plasma pressure and density around 
the x-line drop significantly from the initial sheet value 
(see Supplement A for pressure balance). This fact may 
reduce the impediment that slows the Alfvenic outflows. 
For cases with a higher upstream cr, the initial sheet com- 
ponent is denser and hotter. The system takes a longer 
time to deplete this sheet component and develop rela- 
tivistic inflows, as suggested by the comparison of the 
time-scale between Harris-3 and 5 cases in Table 1. 

In summary, a simple theory based on the Lorentz con- 
traction [12] and the assumption of a universal aspect 
ratio 0.1) of the diffusion region provides an expla- 
nation for the observed relativistic inflows and the en- 
hanced normalized reconnection rate. While the present 
letter was limited to 2D simulations, recent 3D simula- 
tions demonstrate similar relativistic inflows in spite of 
the development of kink instabilities [52] . These results 
may be important for understanding particle acceleration 
[25, 53], the dissipation of strong magnetic field in high- 
energy astrophysical systems, such as the “cr-problem” 
in the Crab Nebula [33], and the destruction of strong 
magnetic field near magnetars [6] and black holes [9]. 

Y.-H. Liu thanks for helpful discussions with S. Zeni- 
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TABLE I: Paramters of Runs 


Harris 

1 

2 

3 

4 

5 

6 

7 

8 

Bg/Bo 

0 

0 

0 

0 

0 

0.2 

0.5 

1 

rib/n ' 0 

1 

0.25 

1 

1 

1 

1 

1 

1 

Tb / me 2 

2.5 

2.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

OJye ! Glee 

0.1 

0.05 

0.1 

0.05 

0.02 

0.05 

0.05 

0.05 

Ox 

4.5 

18.2 

22.2 

88.9 

555.6 

88.9 

88.9 

88.9 

timexccp e a 

200 

100 

250 

500 

1000 

400 

350 

300 

Force-Free 

1 

2 

3 

4 

5 

6 

7 


Bg/Bo 

0 

0 

0 

0 

0 

0 

0 


Tb / me 2 

0.35 

0.36 

0.36 

0.36 

0.36 

0.36 

0.36 


idpe / Gl ce 

1.6 

0.8 

0.4 

0.2 

0.1 

0.05 

0.025 


Ox 

0.1 

0.4 

1.6 

6.6 

26.3 

105.3 

421 



a the time when Vi n reaches a steady high value 


tani, N. Bessho and J. Tenbarge. We are grateful for sup- 
port from NASA through the NPP program and the He- 
liophysics Theory program. Simulations were performed 
at the National Center for Computational Sciences at 
ORNL and with LANL institutional computing. 


[1] H. Ji and W. Daughton, Phys. Plasmas 18 (2011). 

[2] M. Hoshino and Y. Lyubarsky, Space Sci. Rev. 173, 521 

( 2012 ). 

[3] F. V. Coroniti, Astrophys. J. 349, 538 (1990). 

[4] J. Arons, Space Sci. Rev. 173, 341 (2012). 

[5] Y. Lyubarsky and J. G. Kirk, Astrophys. J. 547, 437 

( 2001 ). 

[6] C. Thompson, Mon. Not. R. Astron. Soc. 270, 480 
(1994). 

[7] B. Zhang and H. Yan, Astrophys. J. 726, 90 (2011). 

[8] J. C. McKinney and D. A. Uzdensky, Mon. Not. R. As- 
tron. Soc. 419, 573 (2012). 

[9] K. Beckwith, J. F. Hawley, and J. H. Krolik, Astrophys. 
J. 678, 1180 (2008). 

[10] D. Giannios, Mon. Not. R. Astron. Soc. 408, L46 (2010). 

[11] C. H. Jaroschek, H. Lesch, and R. A. Treumann, Astro- 
phys. J. 605, L9 (2004). 

[12] E. G. Blackman and G. B. Field, Phys. Rev. Lett. 72, 
494 (1994). 

[13] M. Lyutikov and D. Uzdensky, Astrophys. J. 589, 893 
(2003). 

[14] P. A. Sweet, in IAU Symp. in Electromagnetic Phenom- 
ena in Cosmical Physics, ed. B. Lehnet (New York: Cam- 
bridge Univ. Press) (1958), p. 123. 

[15] E. N. Parker, J. Geophys. Res. 62, 509 (1957). 

[16] Y. E. Lyubarsky, MNRAS 358, 113 (2005). 

[17] J. M. Tenbarge, R. D. Hazeltine, and S. M. Mahajan, 
Mon. Not. R. Astron. Soc. 403, 335 (2010). 

[18] S. Zenitani, M. Hesse, and A. Kimas, Astrophys. J. 696, 
1385 (2009). 

[191 L. Comisso and F. A. Asenio, Phys. Rev. Lett. 113, 
045001 (2014). 

[20] M. Melzani, R. Walder, D. Folini, C. Winisdoerfer, and 
J. M. Favre, A & A 570, Alll (2014). 

[21] H. E. Petschek, in Proc. AAS-NASA Symp. Phys. Solar 
Flares (1964), vol. 50 of NASA-SP , pp. 425-439. 


[22] N. Bessho and A. Bhattacharjee, Astrophys. J. 750, 129 

( 2012 ). 

[23] B. Cerutti, G. R. Werner, D. A. Uzdensky, and M. C. 
Begelman, Astrophys. J. Lett. 754, L33 (2012). 

[24] L. Sironi and A. Spitkovsky, Astrophys. J. 783, L21 
(2014). 

[25] F. Guo, H. Li, W. Daughton, and Yi-Hsin Liu, Phys. 
Rev. Lett. 113, 155005 (2014). 

[26] P. A. Sturrock, Astrophys. J. 164, 529 (1971). 

[27] J. Wardle, D. Homan, R. Ojha, and D. H. Roberts, Na- 
ture 395, 457 (1998). 

[28] J. Sakai and T. Kawata, J. Phy. Soc. Japan 49, 747 
(1980). 

[29] A. M. Anile, Relativistic Fluids and Magneto- Fluids 
(Cambridge University Press, New York, 1989), chap. 2, 
p. 34. 

[30] A. Lichnerowicz, Relativistic Hydrodynamics and Mag- 
netohydrodynamics (W. A. Benjamin Inc., New York, 
1967), chap. 4, p. 112. 

[31] Yi-Hsin Liu, W. Daughton, H. Karimabadi, H. Li, and 
S. P. Gary, Phys. Plasmas 21 , 022113 (2014). 

[32] M. Hesse and S. Zenitani, Phys. Plasmas 14, 112102 
(2007). 

[33] J. G. Kirk and O. Skjeraasen, Astrophys. J. 591, 366 
(2003). 

[34] S. Zenitani and H. Hoshino, Astrophys. J. 670, 702 
(2007). 

[35] W. Liu, H. Li, L. Yin, B. J. Albright, K. J. Bowers, and 
E. P. Liang, Phys. Plasmas 18 (2011). 

[36] K. Bowers, B. Albright, L. Yin, W. Daughton, V. Royter- 
shteyn, B. Bergen, and T. Kwan, Journal of Physics: 
Conference Series 180, 012055 (2009). 

[37] Yi-Hsin Liu, W. Daughton, H. Karimabadi, H. Li, and 
V. Roytershteyn, Phys. Rev. Lett. 110, 265004 (2013). 

[38] S. Weinberg, Gravitation and Cosmology: Principle and 
Applications of the central theory of relativity (John Wi- 
ley and Sons, Inc, New York, 1972), chap. 2, p. 51. 

[39] J. L. Synge, The Relativitic Gas (North-Holland Publi- 
cation Company, Amsterdam, 1957), chap. 6, p. 60. 

[40] A. Otto and K. Schindler, Plasma Phys. Control. Fusion 
26, 1525 (1984). 

[41] H. Karimabadi, W. Daughton, and K. B. Quest, Geo- 
phys. Res. Lett. 32, L18801 (2004). 

[42] L. Matteini, S. Landi, M. Velli, and W. H. Matthaeus, 
Astrophys. J. 763, 142 (2013). 

[43] M. Hoshino and K. Higashimori, submitted to J. Geo- 
phys. Res. (2015). 

[44] M. Shigeta, T. Sato, and B. Dasgupta, J. Phys. Soc. 
Japan 54, 3342 (1985). 

[45] S. Wang, L. C. Lee, and C. Q. Wei, Phys. Fluids 31, 1544 
(1988). 

[46] S. Zenitani and M. Hesse, Astrophys. J. 684, 1477 (2008). 

[47] M. Hesse, K. Schindler, J. Birn, and M. Kuznetsova, 
Phys. Plasmas 6, 1781 (1999). 

[48] M. Hesse, M. Kuznetsova, and J. Birn, Phys. Plasmas 
11, 5387 (2004). 

[49] P. Ricci, J. U. Brackbill, W. Daughton, and G. Lapenta, 
Phys. Plasmas 11, 4102 (2004). 

[50] P. L. Pritchett and F. V. Coroniti, J. Geophys. Res. 109, 
A01220 (2004). 

[51] W. Daughton and H. Karimabadi, Phys. Plasmas 14, 
072303 (2007). 

[52] F. Guo, Yi-Hsin Liu, W. Daughton, and H. Li, in prepa- 
ration (2015). 


6 


[53] S. Zenitani and H. Hoshino, Astrophys. J. 562 , L63 

( 2001 ). 


